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Abstract 

Marangoni instabilities in binary mixtures are different from those in pure liquids. In contrast 
to a large amount of experimental work on Marangoni convection in pure liquids, such experiments 
in binary mixtures are not available in the literature, to our knowledge. Using binary mixtures 
of sodium chloride/water, we have systematically investigated the pattern formation for a set of 
substrate temperatures and solute concentrations in an open system. The flow patterns evolve 
with time, driven by surface-tension fluctuations due to evaporation and the Soret effect, while 
the air-liquid interface does not deform. A shadowgraph method is used to follow the pattern 
formation in time. The patterns are mainly composed of polygons and rolls. The mean pattern 
size first decreases slightly, and then gradually increases during the evolution. Evaporation affects 
the pattern formation mainly at the early stage and the local evaporation rate tends to become 
spatially uniform at the film surface. The Soret effect becomes important at the later stage and 
affects the mixture for a large mean solute concentration where the Soret number is significantly 
above zero. The strength of convection increases with the initial solute concentration and the 
substrate temperature. Our findings differ from the theoretical predictions in which evaporation is 
neglected. 
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I. INTRODUCTION 



When a horizontal pure-liquid thin film is sandwiched between a warm solid at the 
bottom and a cool gas at the top, heat transfer can be used to drive fluid flows. At the 
liquid-gas interface, the surface tension is a function of the interfacial temperature, and usu- 
ally decreases with temperature. A small temperature perturbation along the film surface 
may create surface-tension inhomogeneities which in turn trigger fluid instabilities. Flows 
observed in such films are known as Benard-Marangoni convection or surface-tension-driven 
convection. Several recent theoretical and experimental studies are relevant here 
Differing from thin films of pure liquids, the instabilities in binary-mixture films occur 
because the surface tension is a function of both temperature and solute concentration 

g, a 0, a a. 

Given a vertical temperature gradient across the film, the concentration 
gradient may then be imposed by an independent source or generated spontaneously by 
the Soret effect [l(3]. If the thermal Marangoni and solutal Marangoni effects enhance 
each other, a long-wavelength monotonic instability is possible, while if they compete 

nrinriQ 

with each other, an oscillatory instability may occur instead [5|, |6|, 13, la, l9fl. In the past 
decade, theoretical work has led to predictions for the linear and nonlinear stages of 
i—ies in thin filmS of binary natures 0,313,0.3. B y contrast, experiments with 
surface-tension-driven instabilities in binary mixtures are not available in the literature, 
to our knowledge. The lack of such experiments has motivated us to investigate pattern 
formation for binary mixtures in thin films. 



II. GOVERNING EQUATIONS AND PARAMETERS OF THE SYSTEM 



The governing equations of the system, including the Soret effect are: [3, |9| 



V-v = 0, 



d t v + (v • V)u = — p 1 Vp + z/V 2 <7 — g, 
d t T + (v ■ V)T = kV 2 T, 
d t c + (v ■ V)c = D(V 2 c + aV 2 T). 



(1) 

(2) 
(3) 
(4) 
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Here t is time, v, p, p, v and g denote the fluid velocity field, pressure, density, kinematic 
viscosity, and gravitational acceleration, respectively, T and k represent the temperature 
and thermal diffusivity of the fluid, respectively, c is the concentration of the solute, D is 
mass diffusivity of the solute, and a is the Soret coefficient. 

In Eq. ((3D, the heat flux induced by the concentration gradient, known as the Dufour 
effect, is usually very weak in liquids, and thus neglected. At the solid-liquid interface, z = 0, 
v satisfies the no-slip, no-penetration condition v — 0, the temperature is fixed T = T w , and 
the mass flux vanishes, c z + aT z = 0. Here, the wafer temperature T w is constant. At a 



nondeformable liquid-air interface, z = h(t), the boundary conditions are given by [Uj. The 
heat and mass flux balances as well as those for the normal and tangential stresses read, 
respectively, 

k T T z + q(T-T a )+jC = 0, (5) 

- P D{c z + aT z )+jc = 0, (6) 

j = p(w-h t ) = K(T-T a ), (7) 
? 2 

-— = -p + 2pw z . (8) 

Pv 

p,(d z u + Vw) = -a T VT + <j c Vc. (9) 

Here, kx is the thermal conductivity of the liquid, u is the two-dimensional projection of v 
onto the x — y (horizontal) plane, q is the heat transfer coefficient, assuming Newton's law of 
cooling, T a is the constant room temperature, j is the evaporative mass flux, \i is the liquid 
viscosity, p v is the vapor density, and C is the latent heat of evaporation. Also, we have 
assumed that the surface tension a is a linear function of temperature T and concentration 
c; hence, a(T, c) = a — a T (T — T ) + cr c (c — c ). Here, a is the reference surface tension at 
T = T and c = c . (Jt = — (|f )r an d cr c = (f^)c c - AT is the temperature difference across 
the film. 

The dimensionless parameters of the problem are 

P = -L = -,B = \-,M = — , 10 

k T AT k T K aa c 
E= —,K = -r 7 r,X = , H 

respectively, the Prandtl, Lewis, Biot, Marangoni, evaporation, interfacial resistance and 
Soret numbers. The estimated values of these dimensionless parameters at the solute concen- 
tration of 2M/1 (Mole per liter of water) and at the temperature 30°C are P — 7, L^ 1 = 88, 



X = 0.04, E = 2. x 10 -4 , K = 3. x 10~ 3 , B = 6. x 10~ 3 . The Marangoni number depends 
on the instantaneous value of the film thickness h and in our experiments varies between 
2000 and 300. The estimate of the Soret number is based on the experimental data from 



Ref. 



121. 



III. EXPERIMENTAL SYSTEM 



Before carrying out an experiment, the first important issue to be resolved is the selection 
of a resonably non-volatile working liquid, i.e. a binary liquid. We explored a wide range 
of candidate materials. We first tried a mixture of ethanol/water and then abandoned it 
due to the fast evaporation of ethanol. We also tested different combinations of silicone 
oils with noticeable differences in surface tension between them to avoid non-miscibility. 
Only a few silicone oils satisfying the requirements of non-volatility and miscibility are 
commercially available and are affordable. However, for all the possible pairs of oils, either 
one or both are highly volatile at room temperature or they are highly viscous at low 
temperature and volatile when heated. The volatilities are typically much larger than that 
of pure water. As a consequence, we did not choose binary mixtures of silicone oils for our 
experiments. After much trial and error, we chose solutions of sodium chloride (NaCl) with 
pure water as the working binary mixtures. There is an additional advantage in choosing 
the mixture of NaCl/water as a working liquid because its many physicochemical properties 
are available in the literature, including surface tension as a function of both concentration 



and temperature 
concentration 



131 ] - and the variation of the Soret coefficient with temperature and 



12| among others. The evaporation of water is significantly weaker than 



that of other strongly volatile liquids, so that there is a reasonably long time to take 
measurements before the film dries out completely. Finally, the convective patterns are 
readily visible a short time after a thin film is drawn. 



Figured] presents a sketch of the experimental setup of the system. The film was prepared 
on top of a silicon wafer with a circular Plexiglas ring (about 7 cm in diameter) as the lateral 
boundary. The wafer was bonded to an Aluminum disk attached to a thin copper plate. A 
copper tubing coil was soldered to the bottom of the latter. The two ends of the coil were 
connected to a water bath (Neslab RTE-221) using thermally insulated plastic tubing. The 
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temperature of the plate was held constant by water circulation. To further reduce the heat 
loss, the whole copper plate was covered by thermally insulating foam materials, leaving 
only a circular hole to hold the Aluminum disk. The film was directly exposed to the room 
air whose temperature is approximately constant during the course of an experiment. The 
fluctuation in the room air is ±0.5 °C, equivalent to ~ 0.5 x B°C ~ 3mK of the film 
temperature stability. The temperature fluctuation of the wafer was affected by several 
sources. The first source is the thermal fluctuation of the Neslab water bath which does not 
exceed 10 mK. The second source is the temperature fluctuation at the liquid-air interface. 
To further reduce the thermal fluctuation at the silicon wafer, a bridge-controlled heating 



method was applied 



14j. The wafer temperature was measured by a thermistor embedded 



underneath which was measured with an AC bridge. An offset of the bridge balance point 
drove an integro-differential feedback circuit which, in turn, drove a heater attached to the 
inlet of the tubing that carried the cooling water. In such a way, the thermal fluctuation 
of the wafer was further reduced. Depending on the mean temperature of the film, the 
long-time thermal fluctuations of the wafer temperature are less than 30 mK at the mean 
temperature of the film Tf = T w — | A T « 30 °C and less than 15 mK at Tf w 25 °C. The 
mean temperature drop across the film AT is typically less than 1 °C. 



IV. RESULTS AND DISCUSSION 

Figure [2] shows a set of typical shadowgraph flow patterns observed in thin films of 
NaCl/water mixtures with an initial NaCl concentration of AM /I. In this figure, bright and 
dark shades represent downward and upward flows, respectively. In comparison with the 
standard Marangoni convection in pure liquids (MCP), these patterns are very irregular 
and evolve with time. Following the film deposition, an irregular pattern emerges, as seen 
in Fig. 121(a)- This is a transient process with patterns gradually decreasing in size (Fig. 
121(b)). Eventually small-scale patterns cover the entire system, as shown in Fig. 121(c). The 
majority of the patterns are roll-like structures, although there are also localized dots, which 
are rarely observed in MCP. These small-scale patterns persist for a period of time before 
any qualitative changes appear. Then, the characteristic scale of the pattern becomes larger 
and the patterns consist primarily of polygons, as shown in Fig. (21(d-f). If the concentration 
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of NaCl exceeds the maximal solubility (~ 8M/1), salt begins to crystallize out of the 
solution; small islands of salt appear as black dots, as in Fig. E](g-i). An interesting aspect 
of the system evolution concerns the film thickness. Interferometric measurements show 
that except at the meniscus boundary, the film thickness remains horizontally uniform 
during the evolution preceding the formation of salt crystals when deformation of the film 
interface becomes significant. 

To quantify the time variation of the pattern size, we apply fast Fourier transforms (FFT) 
to the original shadowgraph images after the removal of the background intensity variation. 
The images in the Fourier space were averaged over phase angles. From the FFT spectrum, 
we calculate both the mean wavenumber and its standard deviation. The results are then 
converted to a real space pattern size, s(t), as presented in Fig. [3](a). We normalize this 
characteristic pattern size by its instantaneous film thickness, which we refer to hereafter 
as the normalized pattern size. Immediately after the film is initialized, the normalized 
pattern size is approximately 0.8, i.e. the average pattern size is slightly below the value 
of the film thickness. For some time thereafter, the patterns are rolls, as seen in Fig. 
[21(c). The normalized pattern size also retains its initial value for a considerable time, and 
then gradually increases. Due to evaporation, the film thickness decreases monotonically. 
Because the Marangoni number is proportional to the film thickness, it too continuously 
decreases. Finally, the average normalized pattern size reaches s/h ~ 3.0. At this point, 
the corresponding pattern consists mainly of non-equilateral polygons. The irregular shape 



of these polygons leads to large errors in determining the pattern size, Fig. [3j This feature 
is different from the classical MCP where patterns are regular polygons l|, [li]]. This 
difference may be caused by one or a combination of several sources. First, there is the fact 
that the surface conditions in the MCP are usually precisely controlled, as compared to the 
present open system, for which a precise control of the air layer above the film surface is 
much more difficult. Alternatively, we note that the Marangoni number is ever changing, 
and it is impossible to tune the value of the Marangoni number exactly at the threshold value. 



Towards the later stages of the pattern evolution, a noticeable difference in the pattern 
occurs as relatively linear rolls give way to polygonal patterns. At that point, the Marangoni 
number is 539 ± 29, which is quite big compared with the threshold values derived in 
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non-evaporative films [5|, [6fl, and the film thickness was about 0.9 mm. We also carried 
out an experiment with an initial film thickness of 0.87mm, which is close to the value 
mentioned above. At this initial film thickness, the convection starts from polygons and the 
pattern intensity decays with time until the convection pattern vanishes. The magnitude of 
the pattern intensity, I(t), as a function of time is plotted in the inset of Fig. [3](a). I(t) is 
calculated from the average intensity of the image after the removal of the background. A 
large Marangoni number was also observed in an experiment using pure liquid films with 
strong evaporation [ffi| . However, no patterns were observed in pure water films in our 
system. This suggests the solutal effects (thermal solutal and Soret effect) are important. 
This also leads us to the conjecture that a considerable part of the temperature drop 
across the film may actually occur spatially uniformly at the film surface. This is because 
evaporation can create a large temperature drop across the film due to the latent heat loss. 
The rest of the temperature drop across the film is responsible for the fluid instabilities 
at the film surface. Indeed, we estimate that as much as 95% of the temperature drop 
across the fluid may occur at the surface. Note that Fig. [3] does not show the data points 
where the crystallization of NaCl (black dots) takes place at the film surface, as seen in 
Fig. E](g-i). The existence of the solid phase changes the flow pattern dramatically. Solid 
NaCl can grow and merge with a deformable fluid interface. The physics of this stage has 
not yet been well understood and is outside the scope of this work. 

Figure [3](b) shows the time evolution of the normalized pattern size at a different 
temperature, T/ = 25 °C. While the qualitative features are similar to Fig. [3(a), there are 
still some quantitative differences. The normalized pattern size at the early stage is around 
1.0 which is larger than 0.8 at Tf = 30 °C. In the beginning of each run (see e.g. circles and 
diamonds), the normalized pattern size usually decreases for the initial 10 ~ 20mm. This 
feature exists in some of the other runs in figure (a) or (b), except it is more visible here. 
We note that typically, after the film is deposited, patterns generally appear at a slightly 
bigger size; they then gradually decrease in size by developing smaller scale convection 
polygons or rolls such as those in Fig. [2](a-c). 

Figure [3] does not present the normalized pattern-size variation for films with low 
concentrations where the convection flows are weaker and convection decays rapidly. For 
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instance, in Fig. [4](a) the patterns are readily seen at time t = 2.1mm; they are very 
weak if not completely gone in (b) at t = 47.5 min after the film is deposited. Similar 
observations were seen for initial concentrations Co < 0.8 M/l. The decay of the pattern is 
an indication that the local evaporation flux has a tendency to become spatially uniform 
with time. It is reasonable to believe that the evaporation rate, and hence this effect, 
should be independent of concentration. For higher concentrations, however, the patterns 
persist and the intensity tends to be more uniform during the evolution. One possibility is 
that this difference is related to the Soret coefficient, as determined, for instance, from data 
for NaCl/water mixtures by Gaeta et al.[l2| At T = 30 °C, the Soret coefficient is a ~ at 
c = 0.2 M/l and linearly increases as log(c) to a = 1.5 x 10~ 3 /°C at c = 2 M/l. We estimate 
the corresponding Soret number to be x = 0-04 for c = 2 M/l and x = 0-10 for c = AM /I. 
The characteristic time for mass diffusion, about 40 min for a 2 mm thick film, is roughly 
the same as in Fig [U This suggests that the Soret effect becomes important at the later 
stage of convection. In Fig. H](b), the value x ~ 0.002 is too small to drive convection 
after the pattern created by evaporation begins fading away. Thus, in order for the Soret 
effect to play a role, a sufficiently large concentration leading to a higher Soret diffusion 
is necessary since the solutal Marangoni effect is directly proportional to the Soret number X- 



In our experiments we have measured the instantaneous thickness of the film over 
extended periods of time. Using the time history we then evaluate the evaporation rate as a 
function of time. We find that this quantity is nearly constant at a given film temperature 
Tf and a room temperature T a . Considering that the mean concentration increases with 
time, the linear plot for h(t), as seen in the Fig. [5](a), implies a very weak concentration 
dependence on the evaporation rate. The slope of h(t) vs. t increases for a larger film 
temperature Tf or a smaller room temperature T a . The linear decrease of the film thickness 
due to evaporation seems to be robust and almost independent of the NaCl concentration. 
As a comparison, the measurements of the film thickness for evaporating pure-water films 
are also presented in Fig. [5](a). The slope of h(t) for a water film is also a constant, in- 
dicating that this linear behavior is indeed independent of the composition of the liquid film. 

This result is different from the analytical solution obtained by Burelbach et al. [itJ and 
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Oron et al.|ll| for the base state of an evaporating film of a uniform thickness 



h(t) = h ty(K + l) 2 -2Et-K], (12) 

where h is the initial film thickness, K is a constant interfacial resistance to phase change 
and E is the evaporation constant. If Et <C 1, h(t) would be nearly linear. However, as 
Et grows, h(t) would start to bend downwards, especially at small h(t). This trend is not 
seen in the present data. We have experimentally ruled out the possibility that convective 
air motion above the film might be responsible, since the theory assumes that the air layer 
is at rest. In particular, experimental measurements of the evaporation rate for a film at 
a temperature lower than the room air are also presented in the inset of Fig. [5](a) shown 
by the triangle-left symbols; they too show a linear variation. We believe that the reason 
for the discrepancy between theory and experiment is caused by the large fraction of the 
temperature difference that occurs at the liquid-air interface. Because this temperature 
difference remains nearly constant, the evaporation rate also remains constant. Hence, the 
boundary condition used in recent theoretical work differs from that which applies for these 
experiments. 

Keeping other conditions the same, the magnitude of the convection, as determined 
based on the average intensity of the shadowgraph images recorded over the first 15 minutes 
after the film was initialized, depends on the initial concentration c of NaCl and the film's 
mean temperature. At low concentrations, Cq < 0.1 M/l, there are no visible convection 
patterns even for a 2.2 mm thick film. This implies that buoyancy is irrelevant in our 
system if we assume that the thermal expansion of the fluid is not significantly affected 
by adding NaCl to the pure water. When the concentration Co is higher than 0.2 M/l, 
the patterns gradually appear. This observation suggests that thermal Marangoni alone 
cannot drive the fluid instability and that the fluid motion is driven by the combination of 
thermal Marangoni and solutal Marangoni effects. The magnitude increases monotonically 



with concentration for c > 0.2 M/l, At a fixed film temperature of 30 °C , t 
coefficient is negative in the concentration range 0.042 M/l < Co < 0.13 M/l 



le Soret 



How- 



ever, in such a range, the concentration is too low to observe any realistic motions in the film. 
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V. SUMMARY 



To summarize, we have experimentally studied transient Marangoni convection in thin 
binary-mixture liquid films using solutions of NaCl/water in an open system. In the pres- 
ence of evaporation, the patterns, consisting of rolls and polygons, evolve with time. The 
length scale of the pattern increases with time until the convection ceases. The Marangoni 
numbers near onset are much larger than the values predicted in theories developed for 
non-volatile binary liquids. Evaporation is important to the pattern formation at the early 
stage, while the Soret effect is essential at the later stage. For higher concentrations of 
NaCl, the convection becomes stronger. Within the temperature range 25 °C < Tf < 30 °C 
explored in the experiment, convection patterns were observed in the NaCl concentration 
range 0.2 M// < c < 5M/1. Understanding the rich pattern forming dynamics clearly 
requires new theoretical work. 
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Figure 1: Schematic of the experimental setup. Visualization is by standard shadowgraph using an 
expanded laser beam for a light source. The bottom plate of the cell that contains the film consists 
of a highly reflecting silicon wafer that is thermally grounded to an highly conducting copper plate. 
This plate is temperature controlled by means of a regulated water bath. 
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Figure 2: Convection patterns observed at a film temperature Tf = 30 °C and a room temperature 
T a = 20 °C. The initial concentration of NaCl is Co = 4M/I and the initial film thickness is 
ho = 1.8mm. Images (a-i) show the flow pattern at nine different times t. Bright and dark shades 
correspond to downward and upward flows, respectively. The estimated initial Marangoni number 
is M = 2122. 
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Figure 3: The average normalized pattern size s(t)/h(t) as a function of time t. The time axis is 
chosen (thus negative times) so that different data sets have approximately the same film thickness 
h(t) = 0.9 mm at t = Omin. Various symbols represent measurements performed at different initial 
concentrations Co of NaCl, different initial film thickness ho, and different film temperatures Tf. 
The inset shows the normalized pattern intensity I(t)/I(Q) of the data set corresponding to the 
open circles. 




Figure 4: Weakening of the convection intensity at Tf « 30 °C, T a = 22.2 °C, initial concentration 
of NaCl cq = 0AM /liter and initial film thickness ho = 2.2mm. 
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c (M/l) 

Figure 5: (a) The film thickness change as a function of time in the presence of evaporation, (b) 
The convection intensity as a function of the initial NaCl concentration Co and film temperature 7/ 
with the room temperature T a = 20 °C, No convection is observed for cq < 0.2 M/l. 
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